Spectral Distortions at Super-Kamiokande 
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£NJ \ Abstract 

CN ' We examine the effect of the rise in the survival probability of the electron neutrinos 

1 with the decrease in the neutrino energy on the recoil electron spectrum at Super- 

t^- ■ Kamiokande. 



The neutral current (NC) measurements at SNO p] have, conclusively, established the oscil- 
lations of the solar neutrinos and after the evidence for terrestrial antineutrino disappearance 
in a beam of electronic antineutrinos reported by KamLAND [2], all other solutions [31 SJ |5j of 



O 
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SNP can, at best, be just sub-dominant effects. The solar neutrino experiments have, already, 
entered a phase of precision measurements for oscillation parameters. The completeness of 
the LMA solution is being questioned [6] and the scope for some possible sub-dominant tran- 
sitions is being explored vigorously [7J El EJ [Hll EH] ■ The presence of these 'new physics' (NP) 
effects even at a sub-dominant level will affect the present determination of the oscillation 
^ ■ parameters [TT]. As of now, the exact profile of the survival probability of electronic neutri- 

nos over the whole energy spectrum remains unknown as a result of which it is not possible to 
pin point the exact mechanism(s) of neutrino flavor conversion or to exclude the coexistence 
of the other sub-dominant transitions driven by non- standard neutrino- matter interactions/ 
properties. The situation is complicated by the fact that the transition probability in some 
of these non-standard NP scenarios [T2] is energy-independent implying an undistorted solar 
neutrino spectrum which is consistent with the Super-Kamiokande and SNO spectral data. 

In order to study the effect of the rise in the electron survival probability on the recoil 
electron spectrum, we define a quantity S(T) as 



S(T) 



^)lma _ -fir dE H E ) l^ p w + * {1 - P(E)} 



to\ J^dEf(E)^ 



where P(E) is the LMA survival probability [13] and f(E) is the standard 8 B neutrino 
spectrum [TJ]. The quantity S(T) has the interpretation as the probability of an electron 
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being scattered with a recoil kinetic energy T and closely resembles the SK data/SSM ratio 
with the only difference that the SK data/SSM ratio is presented for individual energy 
bins while S(T) as defined above is a continuous function of recoil electron energy. The 
tree level cross-sections for v e e and v x e (x = H,t) scattering have been given in Ref. [15]. 
In fact, S(T) is the LMA expectation for the recoil electron spectrum normalized to the 
standard 8 B neutrino spectrum and has been plotted as a function of T in Fig. 1 along 
with the S(T) for energy-independent asymptotic value of the LMA survival probability for 
comparison. The two probabilities differ considerably from each other and this difference 
could be as large as 10% at 5 MeV. This is enough to highlight the role of spectral distortions 
in discriminating the LMA suppression scenario from the energy independent suppression NP 
scenarios mentioned above [12] . At the current level of precision, both the LMA suppression 
and the energy-independent suppression are consistent with the SK as well as SNO solar 
neutrino data. Consequently, one can not claim a conclusive confirmation of the LMA 
solution from the exclusive study of the high energy region of the solar neutrino spectrum. 
It is, therefore, imperative to study quantitatively the spectral distortions not only to finally 
confirm the LMA solution but also to disentangle the possible NP effects. 

The probability that a recoil electron with true kinetic energy T true will be detected with the 
observed kinetic energy T obs is given by 

r (Ttme , T obs ) = J- exp ( - ( Ttrue I*L. j ^ 

and is called the detector response function. The energy-dependent spread, <tt, is of the 
form 

so that e is the energy spread at 10 MeV. For SK, e = 1.4. To account for the finite energy 
resolution, the differential cross-section must be folded with the detector response function, 
i.e. we must have 




dT t rue ^Tjjy (Ttrue) r (Ttrue, T obs ) (4) 

in Eq. (1) for S(T) to obtain the observed value of 

q/rp , _ /o°° dT true r (T true , T obs ) jgg dEf(E) [^P(E) + P(E)}] 



Jo°° dT true r (T true , T i, s ) IeZZ dEf(E) 



dT tr 



The quantity S(T obs ) has been plotted as a function of T obs in Fig. 2 for the LMA value 
of survival probability as well as the energy independent asymptotic value (EIAV) of the 
survival probability. Also shown are the same two curves depicted in Fig. 1 assuming perfect 
energy resolution. We integrate S(T obs ) over the total recoil electron energy E obs = T obs + m e 
in the bins of 0.5 MeV. This integrated normalized spectrum is denoted by S and has been 
plotted in Fig. 3 alongwith the actual 1496 day SK spectrum normalized to BP04 with 
statistical errors only ^6\. One can calculate the rise in the LMA value SLMA(T obs ) relative 
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to the energy independent asymptotic value SEiAv(T bs) at a particular value of T obs by 
defining 

■n/m \ SLMA\Tobs) ~ SeIAv(T o1)S ) 

R{Tobs) = s • (6) 

iJEIAVy-Lobs) 

It is instructive to see the variation in R(T obs ) because of the energy spread e. Fig. 4 shows 
R(T obs ) for e = (perfect energy resolution) and for e = 1.4 (finite energy resolution). It is 
clear that as a result of finite energy resolution of the detector, Slma{T i, s ) and, therefore, 
R(T obs ) gets enhanced as compared to the corresponding values for perfect energy resolution 
and, again, this increase is more pronounced at the higher energies. Although, the curve for 
R(T obs ) becomes flatter because of the finite energy resolution of the detector but the actual 
value of R(T obs ) becomes larger which is an advantage. One can compare the value of S in a 
low energy bin viz. Sl with that in a high energy bin viz. Sh and find the relative increase 
in S by defining the relative rise-up TZ as 

tz = ^-^. (7) 

<~>H 

It is important to note that TZ defined above is essentially independent of the flux normal- 
ization. Otherwise, it would have been difficult to directly compare S with the experimental 
data which has to be normalized to the SSM 8 B flux which is not known accurately enough. 
Also, the above definition is, almost, independent of 6*13 since the factor of cos 4 #13 cancels 
in the ratio. A non-zero 6*13 will suppress S(T) by the factor of cos 4 6*13 but the 'rise-up' TZ 
will remain practically unchanged (see Fig. 1). Thus, if the SK spectrum is found to differ 
from the LMA spectrum, this conflict cannot be reconciled even in a 3-flavor framework 
since the LMA value of TZ is the same in two/three-flavor framework for small #13 . This 
is contrary to some assertions made recently [17] in literature. Moreover, some part of the 
energy correlated systematic uncertainties will cancel in TZ. 

The variable TZ defined in Eq. (7) can be used as an observable to quantify the turn-up in 
the data/SSM ratio at SK at low energies. It is different from the global observables like the 
moments defined in Ref. [18] in the sense that it directly compares the normalized spectral 
data at two energy ends of the spectrum while the moments are global quantities which are 
not suitable for this purpose. Moreover, the SK solar neutrino data is, already, available in 
small energy bins. For the sake of illustration, we choose the bins E = 5.5 — 6.0MeV and 
E = 13.0 — 13.5MeV and evaluate TZ for these bins to obtain 

TZ = -0.054 ±0.126 (8) 

where both the systematic and statistical errors have been incorporated in the analysis. The 
rise-up is not significantly different from zero for these energy bins. The upper bound on TZ 
is approximately 0.07 (0.20) at la (2a) C.L.. This is to be compared with the corresponding 
LMA value of rise-up 

TZIZa = 0.0871K (9) 

for Am 2 = 7.9 ± 0.3 x 10~ 5 eV 2 and sin 2 9±2 — 0.3^03 where the earth regeneration effects 
have been neglected. However, the earth regeneration effects can be incorporated in a straight 
forward manner which yields 

T^lma = 0-070^ 24 . (10) 
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The rise-up in the total (day+night) SSM normalized rate is 



n LMA = 0.078™. (11) 

This LMA value is compatible with the present experimental value [Eq. (8)] within 1 a 
C.L.. The 'rise-up' TZ in the total (day+night) SSM normalized rate has been plotted in 
Fig. 5 as a function of 812 f° r the central and 1 a upper/lower values of Am\ 2 . Since, the 
rise- up TZ becomes larger for smaller values of 6*12, an upper bound on TZ can be used to 
obtain a lower bound on # 12 Fig. 6 depicts the improvement in the upper bound on 
TZ with increase in statistics and reduction in systematic errors at Super-Kamiokande. The 
confidence levels for the measurement of TZ with the increase in statistics and a projected 
reduction in systematic errors at SK have been given in Fig. 7 from which it is clear that 
with the present level of accuracy at Super-Kamiokande, a rise-up of 20% can be measured at 
l.lcr C.L. Not much improvement results even with the increased statistics with the present 
level of systematic errors. The contribution of hep neutrinos for the bins we have examined 
is very small (about 1.2%) for the SSM value of this flux. However, arbitrarily large values 
of hep flux can substantially affect the higher energy bins with energy greater than 14 MeV. 

The enhancement in S(T) for smaller values of T can be used to further constrain the 
currently allowed neutrino parameter space. We illustrate this point by plotting the constant 
S and TZ curves on the (Amf 2 , 12 ) plane. The quantity S is obtained by integrating S(T) 
over the energy bin of 0.5 MeV centered around T = 5MeV and the quantity TZ is calculated 
from Eq. (7) for the energy bins of 0.5 MeV centered around T = 5MeV, 15MeV. We choose 
these values for illustrative purposes only. The constant S and TZ curves on the (Ara^, 612) 
plane are shown in Fig. 8 for the (Am^, #12) parameter space within the currently allowed 
LMA region. For Am\ 2 and 6 , 12 , we select the 2a ranges of these quantities and plot the 
constant S curves for S = 0.44, 0.46, 0.48 (curves with negative slope from left to right, 
respectively) and, also, plot the constant TZ curves for TZ = 0.12,0.10,0.08 (curves with 
positive slope from left to right, respectively). It is evident that there is an increase in S 
and decrease in TZ with increasing 12 . Therefore, an accurate measurement of S and TZ at 
5MeV will further constrain the LMA allowed 612- For instance, if S is found to be smaller 
than 0.46 (i.e. S < 0.46) and TZ is found to be less than 12% (i.e. TZ < 12%), 8 12 will, 
approximately, be within the range 31° < 812 < 34°. In conclusion, the prospects for the 
observation of spectral distortions at SK-III, in case it comes up, appear to be bright [19J. 
The failure to observe spectral distortions at SK-III would signal new physics beyond LMA. 
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Figure 1: S(T) versus T for 13 = 0° (left panel) and 13 = 12.3° (right panel). 
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Figure 2: Probability folded with detector response function and neutrino cross-sections as 
the function of recoil electron kinetic energy. 
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Figure 3: The recoil electron spectrum normalized to SSM in the bins of 0.5 MeV (S). The 
SK data with statistical errors is also shown for comparison. 
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Figure 4: R(E obs ) versus E obs . 
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Figure 5: The rise-up 1Z in the total (day+ night) SSM normalized rate versus 612 for the 
central and 1 o upper /lower values of Am 2 . The solar mixing angle #12 has been varied over 
its 1 a range. 




Figure 6: Improvements in the upper bound on TZ with increased statistics with present (left 
panel) and half of the present (right panel) systematic errors. The solid, dashed and dotted 
lines are the upper bounds at 1, 2 and 3 a C.L. 
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Figure 7: (a) Confidence levels for the measurement of 1Z with increase in statistics for 
present systematic errors (left panel), (b) Confidence levels for the measurement of TZ with 
increase in statistics for half of the present systematic errors (right panel). The solid, dashed 
and dotted lines are for TZ = 0.10, 0.15 and 0.20. 
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Figure 8: Constant S and 1Z plots. 
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